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基于 人 工 鱼 群 算法 的 分 数 阶 P 上 控制 器 参数 整定 
张 学 典 ， 王 富 庆 ， 秦 晓 全 


(上 海 理工 大 学 光电 信息 与 计算 机 工程 学 院 , 上 海 200093) 


摘 要 : 针对 分 数 阶 PT 控制 器 的 参数 整定 ， 提 出 图 像 法 和 人 工 鱼 群 优化 算法 相 结合 的 方法 来 对 分 数 阶 PT 控制 器 进 
行 参数 整定 。 分 别 以 一 阶 系统 和 二 阶 系统 来 模拟 典型 的 速度 0 统 ， 以 此 模型 为 被 控 对 象 进行 分 数 阶 PP 控制 器 的 
设计 ， 首 先 在 频 域内 ， 根 据 系 统 的 相对 稳定 性 和 增益 变化 的 鲁 棒 性 等 条 件 ， 推 导出 方程 式 ; 然后 根据 图 像 法 解 出 分 数 
阶 pT* 控制 器 的 参数 ， 以 解 出 的 参数 为 中 心 位 置 ， Sap 进而 用 人 工 鱼 群 算法 对 其 周围 进行 寻 优 。 最 后 进 
行 仿真 研究 ， 通 过 仿真 可 以 看 出 ， 通 过 人 工 鱼 群 算法 寻 优 得 到 的 控制 器 ， 比 单纯 用 图 像 法 得 到 分 数 阶 PI 控制 器 能 使 
系统 具有 更 好 的 动态 响应 特性 ， 并 且 满 足 增益 变化 鲁 棒 性 的 条 件 。 
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Parameter tuning of fractional order PI* controller based on artificial fish school algorithm 
Zhang Xuedian, Wang Fuyan, Qin Xiaofei- 


(School of Optical-Electrical & Computer Engineering, University of Shanghai for Science & Technology, Shanghai 200093, 
China) 


Abstract: For parameter tuning of fractional order py* controller, this paper proposed the artificial fish swarm algorithm 


combined with the image method to tune the parameter. The paper used the typical first-order and second-order system to 


represent the typical speed servo system and used this model to set the parameters of the fractional order py* controller for 


和 的- the controlled object. Firstly, in the frequency domain, according to the relative stability of the System and the robustness of 


the gain, the paper derived the equation. After that, using the image method , it can solve the parameters of the fractional 


order pI* controller Around the parameter, according to artificial fish swarm algorithm, it can optimize these parameters. 
Finally, simulation results Show that the controller optimized have better dynamic characteristic than controller obtained by 


the image method and meet the conditions of robust gains. 
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可 调 参 数 4 和 /， 使 控制 器 具有 更 好 的 灵活 性 四。 
传统 的 整数 阶 PID 的 整定 方法 已 经 得 到 了 很 多 深入 的 丰 
最 近 分 数 阶 微 积分 逐渐 引起 国内 外 学 者 的 关注 ， 尤 其 在 系 。” 究 63， 在 此 基础 上 很 多 学 者 提出 了 分 数 阶 PPpDr 控制 器 整定 

统 理论 领域 和 控制 领域 ， 因 为 可 以 用 分 数 阶 搭建 更 准确 的 数学 方法 。 胡 海 波 等 在 传统 的 粒子 群 的 基础 上 ， 使 细菌 员 食 算法 与 

模型 ， 增 强 系统 的 控制 性 能 023。 分 数 阶 微 积 分 的 概念 已 经 应 ”粒子 群 算法 相 结合 对 分 数 阶 PID 控制 器 进行 设计 ， 提 高 了 传 

用 到 很 多 工程 领域 ， 从 而 产生 了 很 多 分 数 阶 控制 器 。 ” 统 粒 子 群 算法 的 计算 效率 ， 同 时 得 到 的 分 数 阶 PID 控制 器 ， 

Oustaloup 以 控制 动态 系统 的 思想 ， 开 发 了 CRONE 控制 器 ， ” 比 整 数 阶 控制 器 具有 更 好 的 动态 特性 外。 秦 君 琴 等 人 四 把 分 数 

且 证 明了 CRONE 控制 器 比 传统 的 整数 阶 PID 控制 器 具有 较 。 阶 PID 控制 器 应 用 到 大 棚 室 内 的 温度 调节 ， 使 其 产生 的 效果 

大 的 优势 四 。Matignon 等 人 针对 分 数 阶 系统 的 稳定 性 、 可 控 性 与 传统 的 整数 阶 PID 和 模糊 PID 控制 器 进行 对 比 ， 最 后 发 现 

和 可 观 性 进行 了 研究 。 随 后 ，Podlubny 在 研究 分 数 阶 控制 系 ”分 数 阶 控制 器 提高 了 系统 的 自 适应 性 和 和 鲁 棒 性 。 同 时 ， 

统 基础 上 在 1994， 提 出 了 分 数 阶 PPDy 控制 器 。 相 比 于 传统 。 Mohanty 等 人 no 针对 独立 微 电 网 的 无 功 补偿 ， 利 

的 整数 阶 PID, 上 有 具有 更 好 的 鲁 棒 性 和 准确 性 ， 并 且 增 加 了 两 个 ICA(imperialist competitive algorithm) 进 行 分 数 阶 PID 控制 器 
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的 设计 ， 最 后 通过 改变 输入 参数 和 负载 条 件 ， 来 验证 分 数 阶 控 
制 器 的 优势 。Jauregui 等 人 0 用 改进 的 粒子 群 算法 对 分 数 阶 
PPDa 控制 器 进行 了 整定 。Altintas 等 人 023 通 过 用 遗传 算法 对 
磁悬浮 系统 进行 控制 器 的 参数 的 优化 ， 从 而 获得 分 数 阶 PP 
空 制 器 (FOPID) 和 整数 阶 PID 控制 器 (IOPID)， 通 过 分 析 验 证 
了 FOPID 具有 更 好 的 灵活 性 ， 使 系统 具有 很 好 的 动态 响应 特 


ChinaXiv 合 作 期 刊 


张 学 典 ， 等 : 基于 人 工 鱼 群 算法 的 分 数 阶 控制 器 参数 整定 


根据 控制 系统 增益 的 鲁 棒 性 条 件 ,以 增益 鲁 棒 性 为 ? 
作为 控制 器 的 设计 的 约束 条 件 ， 增 益 鲁 棒 性 准则 如 式 (6) 所 示 。 

dLArglCdwPUo) 
cdw 0 


增益 鲁 棒 性 准则 规定 在 增益 穿越 频率 w 处 ， 系 统 相 频 特 


住 则 [17] 


=0 (6) 


[ 


性 曲线 导数 为 0，， 即 系统 的 相 频 曲线 是 水 平 的 ， 当 系统 增益 发 


性 。Singh03l 设 计 了 基于 蚁 群 算法 分 数 阶 模糊 pppr 控制 器 。 
Jain 等 人 04 在 直流 电机 速度 控制 环节 ， 用 粒子 群 算法 设计 了 分 
数 阶 pppr 控制 器 。 

本 文 主要 是 在 文献 [15] 运 用 的 图 像 的 基础 上 ， 运 用 人 工 鱼 
群 算法 对 分 数 阶 pm 控制 器 三 个 参数 进行 优化 ， pm 控制 器 的 
型 公式 ， 如 式 (D 所 示 。 


C(s)=K,(+ (1) 
Ss 


其 中 :,,k 是 传统 的 整数 阶 可 调节 的 参数 ，4 为 积分 阶 次 。 

人 工 鱼 群 算法 是 有 李 晓 厦 2002 年 提出 的 na, 他 是 一 种 基于 
动物 行为 的 智能 优化 方法 ， 对 寻 优 的 空间 和 性 质 没有 特殊 要 求 ， 
这 种 搜索 方法 通过 比较 目标 函数 值 来 得 到 结果 ， 不 需要 有 导数 
信息 ， 因 而 具有 较 好 的 全 局 寻 优 能 力 ， 可 以 并 行 搜索 ， 搜 索 速 
度 快 等 优点 。 由 于 人 工 鱼 群 算法 对 初 值 不 敏感 ， 避 免 寻 优 过 程 
中 发 生 局 部 最 优 的 情况 。 

分 别 用 图 像 法 和 人 工 鱼 群 算法 得 到 分 数 阶 Pr 控制 器 的 参 
数 ， 然 后 通过 仿真 ， 对 比 两 种 方法 的 仿真 结果 ， 可 以 发 现 用 人 
工 鱼 群 优化 的 结果 相 比 图 像 法 ， 能 使 系统 具有 更 好 的 动态 响应 
特性 ， 最 后 使 此 优化 结果 与 其 周围 的 参数 进行 比较 ， 仿 真 图 像 
表明 优化 的 结果 能 使 系统 达到 最 佳 的 状态 。 


1 ”PF 图 像 法 的 实现 


个 可 调 


分 数 阶 Pr 控制 器 相对 于 整数 阶 PI 控制 器 多 了 参 
数 4， 能 在 更 大 范围 内 调节 控制 器 的 性 能 ， 所 以 采用 分 数 阶 
pr 控制 器 的 系统 能 达到 更 好 的 控制 性 能 。 分 数 阶 Pr 控制 器 
设计 准则 要 考虑 稳定 性 和 增益 鲁 棒 性 的 要 求 。 以 典型 的 一 阶 系 
统 P(s)， 如 式 (2) 所 示 ， 作 为 被 控 对 象 ， 用 图 像 法 设计 其 分 数 


阶 PE 控制 器 。 


PO 2 

基于 分 数 阶 Pr 控制 器 典型 的 一 阶 惯性 系统 的 开 环 传递 函 
G(s) = C(s)P(s) (3) 

其 中 : C(s) 是 典型 的 分 数 阶 PT 控制 器 如 式 (1)。 根 据 系 统 的 相 


频 特性 和 幅 频 特性 ， 首 先 指定 系统 的 截止 频率 w ， 在 截止 频 
率 处 幅 值 满足 式 (4)。 


ICcGw)PGw)|=1 (4) 
为 了 保证 系统 的 稳定 性 ， 指 定 其 相位 裕 度 om, ， 可 列 出 下 
列 方程 式 : 


Arg[CG@®)PI@)]=-x+9, G5) 


生 小 的 扰动 时 ， 系 统 的 相位 裕 度 基本 保持 不 变 ， 在 阶 跃 响应 
线 上 通过 观察 他 的 超 调 量 ， 看 是 否 满足 增益 鲁 棒 性 的 条 件 ， 此 
时 超 调 量 应 该 基本 保持 不 变 09。 
对 式 (4)~(6) 进 行 化 简 ， 可 以 得 到 下 列 公式 ; 


k, d+ ko, ?cos(A5/)) + @, sin( A457))? 
|GGowj = 一 7% % =1 


> 


1+(@T) 
(7) 
| | 加 天 Ma sin(22) 
dw jo 0. +2k08 cos(h)+k, (8) 
I 
1+(T@.) 
。 1, kw sin(s) 了 
4 = Ut T)= 
Bien 机 大 OO ent WW (9) 
-NX+9, 
联 立 (7),(8),(9) 三 个 方程 组 解 得 (10) 和 (11): 
| -tan[tan  (@. T) + 9,] da0) 
四 一 Sin(2 芒 ) 十 四 CoOS(22) tan[tan (mw T)+9,] 
4o 人 +B +TA=0 (11) 
由 式 (11) 得 : 
BYB -AA i 
24w 
其 中 
B=2A@ 2 cos(02) 一 sin(22) 
了 
1+CToo) 


给 定 截 止 频率 、 相 角 裕 度 、 时 间 常数 w =10rad/s， w = 
70 ，7=0.6s。 分 别 画 出 式 (10)(12) 的 图 形 ， 关于 4 的 函数 
像 ， 如 图 1 所 示 。 


| 


80 7 
一 一 (10) 
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图 1 图 像 法 获得 的 曲线 


12 14 16 18 2 


红线 是 式 (12) 的 波形 图 ， 蓝 线 是 式 (10) 的 波 


DS 
砚 


， 则 它 
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们 的 交点 为 所 求 的 参数 ， 则 (4, ) 的 值 为 (0. 


72,14.75)， 代 入 


式 (7) 可 以 得 到 k, =2.868 进而 得 到 分 数 阶 PT 控制 器 : 


CG) =2.868( + ) 


G(s) 的 闭环 传递 函数 ， 对 其 P 
则 可 以 得 到 系统 的 阶 跃 响 应 


到 


砚 
US 


攻 。 如 图 


2 所 示 


Step response 


式 (3) 可 以 求 出 ， 系 统 的 开 环 传递 函数 G(s) ， 然 后 求 出 
环 传递 换 函 数 进 行 阶 跃 仿 真 ， 


(13) 


o 


0 0.5 1 ee 2 2.5 3 
图 2 图 像 法 得 到 的 系统 阶 跃 响应 
可 以 从 阶 跃 响应 曲线 中 得 出 系统 的 性 能 指标 。 一 个 是 上 升 
时 间 7 =0.13， 调 节 时 间 + =0.82， 超 调 量 o% =19%。 阶 跃 响 
应 的 上 升 时 间 是 系统 从 0 起 第 一 次 上 升 到 稳 态 值 的 时 间 ， 反 映 
了 系统 的 响应 的 快速 性 ， 调 节 时 间 是 系统 到 达 稳 态 值 附近 的 最 


| 


j 节 的 快慢 。 超 调 量 是 输出 量 超 出 稳 态 值 的 最 


短 时 间 ， 一 般 是 稳 态 值 2% 以 内 ， 这 里 取 0.998， 它 衡量 系统 


大 偏 移 量 与 稳 态 


值 之 比 ， 如 式 (14) 所 示 ， 它 反映 了 系统 的 相对 稳定 性 。 图 像 法 
是 通过 指定 截止 频率 w 来 求解 41 和 ， 进 而 求 得 k, ， 不 能 随 


意 给 定 截止 频率 w ， 从 而 得 不 到 系统 的 最 优点 ， 不 能 


工作 在 最 优 状态 。 因 此 使 用 智能 优化 算法 和 系 


使 系统 
统 的 增益 鲁 棒 性 


结合 ， 对 (wo ,42) 参数 进行 优化 。 系 统 的 截止 频率 w 不 会 特定 


给 出 ， 使 得 w 和 多 自由 
值 ， 从 而 得 到 最 优 的 组 合 。 


组 合 ， 用 人 工 鱼 群 算法 筛选 出 最 优 的 
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b) 对 准 规则 ,尽量 与 临近 伙伴 的 平均 方向 一 致 。 


c) 内 聚 规则 ， 

追尾 行为 是 
鱼 群 则 会 跟着 向 
水 域 中 食物 
状态 值 为 待 优化 
为 ， 通 过 这 些 方 
据 食 物 的 浓度 随 


尽量 朝 临 近 伙 伴 的 中 心 移动 。 


指 当 鱼 群 中 某 个 鱼 发 现 食物 多 的 水 域 ， 则 其 他 


该 水 域 游 动 。 
的 浓度 是 人 工 鱼 群 算法 的 
的 变量 。 
式 ， 可 以 得 到 最 优 的 解 。 在 砚 


保证 了 鱼 群 的 多 样 性 。 


聚 群 行 


时 


人 工 鱼 与 周围 环 


| 


境 相互 作用 的 过 程 ， 这 两 种 行为 保 订 


我 学 习 的 过 程 ， 
文 


间 不 太 拥挤 ， 可 
而 寻找 到 最 有 的 
鱼 群 通过 自身 的 
的 平衡 。 
人 工 鱼 群 的 


最 终 寻 


标 函 数 ， 每 个 鱼 的 
它 包 括 疯 食 行为 、 聚 群 行为 和 追尾 行 
食 行为 中 鱼 群 根 
机 移动 ， 它 属于 极 值 优化 的 过 程 ， 同 时 也 是 自 
为 和 追尾 行为 时 
E 人 工 鱼 之 


以 保证 人 工 鱼 总 体 朝 着 食物 多 的 地 方 游 动 ， 从 


伙伴 。 人 工 鱼 是 一 种 群体 智能 优化 算法 


。 人 工 


信息 和 周围 环境 来 调节 移动 的 方向 ， 实 现 种 群 


区 到 食物 浓度 最 大 的 人 工 鱼 的 位 置 


如 果 没 有 进步 在 


三 种 模式 可 以 互相 嵌 套 ， 可 以 先进 行 追尾 行 


为 
要 


进行 名 


则 执行 聚 群 行为 
模式 最 终 的 结果 
则 鱼 群 就 会 向 最 
程 中 就 有 下 面 的 
a) 并 行 性 ,多 


， 如 果 仍 然 没 有 进步 则 就 执行 
就 是 向 系统 最 优 的 方向 前 进 ， 
大 值 方 向 移动 。 所 以 人 工 鱼 群 
特征 : 
个 AF 并 行 的 进行 搜索 。 


随机 行为 


b) 简 单 性 ,算法 


中 仅 使 用 了 目标 问题 的 函数 值 。 


c) 全 局 性 ,算法 具有 很 强 的 跳出 


食 行 为 中 具有 随 


局 部 极 值 的 能 力 ( 医 


食 行为 ， 如 果 还 没有 朝 指定 方向 前 进 ， 
。 三 种 
如 果 是 求 最 大 值 
算法 在 寻 优 的 过 


为 在 贡 


机 性 )。 


d) 快 速 性 ,算法 中 虽然 有 一 定 的 随机 因素 ， 但 总 体 是 步 步 


向 最 优 搜索 。 


e) 跟 随 性 , 随 着 工作 情况 或 其 他 


的 漂移 ， 具 有 快 


图 
测 


速 跟踪 变化 的 能 


2.2 人工 鱼 群 模型 的 构造 
人 工 鱼 群 算法 采用 自 下 而 上 的 设计 方法 : X = (xx …,X ) ， 


其 中 二 G=12…， 


m) 为 寻 优 的 变量 ， 人 


示 为 FC=F(9)， 
为 d; | XX | 


gen 为 迭代 次 数 ， 
当 fishnum 


办 代 次 数 就 会 增 


的 变更 造成 的 极 值 点 


笛 当 前 的 食物 浓度 表 
其 中 FC 是 目标 函数 ， 人 工 鱼 个 体 之 间 的 距离 


; visual 表示 人 工 鱼 的 感知 距离 ; step 表示 人 
工 鱼 移动 的 最 大 步 长 ; 6 为 拥挤 度 因 子 ，fishnum 为 鱼 的 数量 ， 


try_number 为 尝试 次 数 。 
越 大 有 利于 发 生 局 部 最 优 ， 收 敛 速度 很 
加 。 当 求 最 大 值 问 题 时 ， 取 0<6<1, 当 


时 ， 它 的 允许 
但 准 


a% = Cu 一人 (14) 
Cu 
2 “人 工 鱼 群 算法 
2.1 人 工 鱼 群 的 特点 
人 工 鱼 群 算法 是 一 种 基于 动物 行为 的 群体 智能 优化 算法 ， 
这 种 算法 来 源 于 鱼 的 疯 食 行为 。 由 于 鱼 的 自主 性 ， 鱼 往往 能 自 
行 或 尾随 其 他 鱼 ， 找 到 食物 源 最 多 的 地 方 。 因 此 鱼 数目 多 的 地 
方 就 是 本 水 域 食 物 多 的 地 方 。 根 据 这 些 特点 ， 来 构造 人 工 鱼 群 
算法 ， 通 过 人 工 鱼 来 模拟 鱼 群 的 砚 食 、 聚 群 、 追 尾 和 随机 行为 ， 
从 而 实现 问题 的 寻 优 。 下 面 列举 了 鱼 的 三 种 典型 的 行为 : 
况 食 行为 :， 当 鱼 在 水 域 中 自由 移动 时 ， 当 鱼 发 现 食物 多 的 


区 域 ， 鱼 会 自主 向 该 区 域 游 动 。 
聚 群 行为 : 为 了 保证 生存 和 躲避 伤害 ， 鱼 
。 聚 群 行为 遵守 三 条 规则 ; 


会 自然 地 聚集 成 


a) 分 隔 规 则 ,尽量 避免 与 临近 伙伴 过 于 拥挤 。 


允许 拥挤 度 也 会 
的 准确 度 也 会 降 
行为 就 越 显 著 ， 


确 度 会 降低 。 


挤 度 就 越 小 ， 在 整个 范围 内 
如 果 求 最 小 值得 问题 ， 


取 6>1， 当 


快 ， 但 


6 越 大 


kt 有 更 好 的 收敛 性 ， 
0 越 小 ? 


变 小 ， 在 整个 区 域 中 有 很 好 的 收敛 性 ， 同 时 它 
低 。 当 视野 范围 越 小 ， 则 融 食 行为 和 随机 移动 
当 视 野 范围 visual 越 大 时 ， 妃 尾行 为 就 越 显著 。 


step 越 大 ， 有 利 
震荡 ， 影 响 后 期 


寻 优 ， 但 在 整个 


于 早期 的 寻 优 ， 但 在 后 


期 收敛 过 程 中 容易 产生 


收敛 的 速度 。step 越 小 ， 


利于 在 部 分 


x 域 的 


x 域 的 寻 优 结果 变 慢 。 


]) 竞 食 行为 
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设 入 工 鱼 当 前 状态 为 Xx, ， 在 其 感知 的 范围 


内 随机 选择 一 


个 状态 X, ， 如 果 求 极 大 值 问题 时 ， 当 FC, > FC, ， 则 可 以 向 


该 方向 


Xi， 后 了 


些 


断 是 否 满足 前 进 的 


前 进一步 ， 如 式 (15) 所 示 。 反 之 ， 则 再 重新 选择 状态 


要 求 ， 如 果 反 复 尝 试 try_number 


次 后 ， 如 果 不 满足 要 求 则 随机 移动 一 步 ， 如 式 (16) 所 示 。 


Step*rand() 取 [0,step] 之 间 的 任意 的 值 ， 


X,(next) = X, a 


X,(next) = 


,一 X， 
— * step* rand() 


个 鱼 的 位 置 移动 。 


2) 聚 群 行为 


人 工 鱼 的 当前 状态 为 XX ， 探 索 当前 
_，FC 代表 伙伴 中 心 位 置 事物 的 


的 伙伴 数 


浓度 。 当 FC./n, > 5FC, ， 则 表明 伙伴 中 心 有 较 多 的 


目 及 中 心 位 置 X 


| (15) 

xX, 一 Xi 

X,+ step * rand() (16) 
一 个 阐 刷 二 
| 和 ij 一 Xi | 


领域 内 ( d; < visual ) 


食物 


且 


不 太 拥 挤 ， 则 朝 伙伴 的 中 心 位 置 方向 前 进一步 ， 则 执行 式 (17)， 
否则 执行 更 食 行为 。 
pe X.—X, 米 米 
X,(next) = X,+ Xx step * rand() (17) 
3) 追 尾行 为 
人 工 鱼 当前 状态 为 X,， 探索 当前 领域 d; < visual ) 有 最 
大 实 食 物 浓度 FC， 人工 鱼 的 位 置 X%，。。 如 果 
FC 1n >6FC,， 表 明 伙伴 Xx 的 状态 具有 较 高 的 食物 浓度 


且 不 太 拥 挤 ， 则 朝 X 方向 前 进 


X,(next) = X,+ 


4 公告 板 


公告 板 是 上 


处 位 


上 的 值 进 行 比较 ， 检 验 自 


于 公告 板 ; 


分 (ITAE 


和 选择 性 的 控制 系统 性 能 评价 指标 09]， 


上。 -Xl 


步 ， 则 执行 式 (18), 否 则 


一 了 X， 
ee ~ * step* rand() 


(18) 


来 记录 最 优 人 工 鱼 个 体 的 状态 和 和 该 人 工 鱼 所 
的 食物 浓度 ， 每 次 优化 完毕 后 ， 使 优化 的 结果 与 公告 板 
身 状 态 与 公告 板 的 状态 ， 如 果 自 身 优 


大 态 ， 则 替换 公告 板 上 的 最 优 状 态 。 当 达到 迭代 次 数 
后 ， 公 告 板 上 显示 的 值 为 最 优 值 。 
2.3 人工 鱼 群 的 设计 及 其 流程 


在 对 分 数 阶 PT 控制 器 进行 优化 过 程 中 ， 
标 函 数 ， 即 食物 浓度 ， 本 次 优化 实验 中 


需要 构造 一 个 
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鱼 的 位 置 X 可 以 表示 为 X(wo,,4), 人 工 鱼 之 间 的 距离 ， 如 式 
(20) 所 示 。 


d; = VY(@®,; 一 四) +(Mi -4) (20) 


i 


yd 


通过 限制 w 和 4 的 范围 ， 缩 小 寻 优 的 时 间 ， 增 加 准确 度 。 
以 图 像 法 得 到 的 和 入 作为 中 心 位 置 。 设 定 它们 的 取 值 范围 ， 
如 式 (21) 所 示 。 


Gd-ano <w.<(I+a)o. QD 
d=-a4 <4<(+an1 
其 中 0<a<1。 


设计 的 算法 流程 图 如 图 


3 所 示 。 


用 图 像 法 求 出 寻 优 范围 的 中 心 
值 〈 确 定 寻 优 的 范围 


! 


初始 化 鱼 群 ， 确 定 各 个 变量 的 
值 。 


J 


在 给 定 的 范围 内 随机 生成 各 个 

鱼 的 状态 X， 同 时 计算 出 各 个 

状态 对 应 的 目标 函数 值 Fc。 并 
给 公告 板 赋值 。 


用 鱼 群 中 的 三 种 行为 进行 筛 
选 。 得 到 三 者 中 最 优 的 个 体 X 
和 其 对 应 位 置 的 食物 浓度 FC. 


与 公告 板 上 的 值 进 行 比较 。 是 
否 大 于 公告 板 上 的 值 。 


~ 


< 


保留 公告 板 上 
的 值 X 和 FC 


替换 公告 板 


上 的 FC 和 Xx. 有 一 | 


(E 千 极 上 的 人 为 最 
优 值 


\ 


6/ 


图 3 算法 流程 


人 工 鱼 群 算法 求 分 数 阶 PT 控制 器 详细 的 寻 优 步骤 如 下 : 

a) 初 始 化 鱼 群 ， 确 定 个 变量 的 值 ， 取 人 工 鱼 数量 
fishnum=100， 最 大 和 连 代 次 数 M4AXGEN=50， 设 定 人 工 鱼 的 可 
视 域 visual=1， 步 长 step=0.01， 拥 挤 度 因子 56=0.618， 尝 试 次 
数 try_number=100。 在 式 (21) 的 范围 内 随机 生成 fishnum 个 人 
工 鱼 的 状态 X,(o,,?) 。 通 过 式 (7) (10) 计 算 各 个 系数 ， 最 后 计 


以 时 间 误差 绝对 值 积 


) 作为 目标 函数 ， 


使 函数 JQTAE)=[ +leOld 的 值 最 小 ， 在 本 文中 通过 求 


已 是 一 种 具有 很 好 的 工程 实用 性 


在 进行 优化 过 程 中 


JITAE) 倒数 作为 目标 函数 如 式 (19)， 即 用 人 工 鱼 群 算法 求 出 
它 的 最 大 值 。 

es (19) 

确定 了 目标 函数 以 后 ， 还 要 确定 他 的 搜寻 范围 ， 这 样 可 以 


日 


结果 的 精度 ， 降 低 搜寻 的 时 间 。 在 本 文中 ， 每 个 人 工 


提高 搜寻 结果 


算出 FC=1/JGTAE) 。 此 时 0< w, <20,0< 4 <1.42。 同 时 给 公告 
板 赋 初 值 。 


、 追 尾 和 聚 群 行为 ， 搜 寻 具 有 最 大 食物 浓 
EE， 得 到 新 的 个 体 xX, 和 FC。 

c) 对 通过 比较 三 种 行为 得 到 的 个 体 X 和 FC， 得 到 三 者 中 
最 优 的 个 体 ， 再 与 公告 板 上 的 值 进行 比较 ， 若 该 位 置 的 食物 浓 
度 FC 大 于 公告 板 ， 则 替换 公告 板 上 的 人 工 鱼 位 置 X 和 食物 
浓度 FC。 

d) 当 达到 最 高 迭代 次 数 MAXGEN， 算 法 结束 ， 此 时 的 公 
告 板 上 的 值 X(o,,4) 要 求 得 最 优 值 。 如 果 gen 没 达 到 最 大 值 ， 
则 gen=gent+1， 转 b)。 


法 
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3 ”仿真 检验 


典型 的 速度 环 控制 模型 如 图 4 所 示 。 


图 4 典型 的 速度 环 控 制 模型 


人 工 鱼 群 算法 迭代 过 程 如 图 5 所 示 。 


鱼 群 迭代 过 程 


5 
加 


口 
和 


日 标 函数 们 


和 代步 


群 迭 代 过 程 


图 5 人 


通过 图 5 可 以 看 出 ， 在 和 欠 代 20 次 左右 时 ， 目 标 函 数 基本 
达到 了 稳定 状态 。 和 迭代 20 次 后 已 经 确定 了 最 优 值 ， 说 明 寻 优 
的 过 程 还 是 很 快 的 。 
P(s) 为 典型 的 一 阶 速度 模型 ，C(s) 为 分 数 阶 PT 控制 器 ， 
通过 人 工 鱼 群 算法 优化 后 得 到 w=16.43， 4=0.58。 则 k, =29， 
代入 式 (7) 可 以 得 到 k, =3.487。 分 数 阶 PT 控制 器 为 


CG)=3.4874+ 熏 QD 


代入 式 (3) ， 同 样 得 到 系统 的 开 环 传递 函数 G,(s) 。 为 了 
验证 优化 后 系统 的 动态 特性 得 到 改变 ， 通 过 阶 跃 响应 曲线 对 比 
两 种 方法 的 动态 特性 。 同 时 在 w=16.43，44.=0.58 周围 选取 四 
个 代表 性 的 值 ， 形 成 四 个 分 数 阶 Pr 控制 器 ， 与 用 人 工 鱼 群 算 
法 的 得 到 的 控制 器 进行 阶 跃 响应 的 比较 。 为 了 验证 增益 鲁 棒 性 ， 
分 别 取 优 化 后 k, 的 100%，110%，90%， 对 比 不 同 k, 的 情况 
下 系统 的 阶 跃 响 应 。 

首先 对 用 两 种 方法 得 到 的 控制 器 进行 比较 ， 分 别 得 出 
G(s)、 G,(s) 对 应 的 闭环 传递 函数 的 阶 跃 响应 曲线 ， 如 图 6 所 


不 。 


EA 


Step response 


G2 


0 0.5 1 1.5 2 2.5 3 
Time [s] 


图 6 图 像 法 与 人 工 鱼 群 算法 比较 
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通过 图 6， 把 两 个 系统 的 阶 跃 响应 的 性 能 指标 以 表格 的 
式 进行 对 比分 析 ， 如 表 1 所 示 。 


SS 


表 1 两 种 方法 指标 对 比 


调节 时 间 上升 时 间 超 调 
Gl 0.82 0.13 19% 
G2 0.36 0.07 16.6% 


通过 表 1 可 以 看 出 ， 优 化 得 到 的 分 数 阶 PT 控制 器 ， 调 节 
时 间 ， 上 升 时 间 和 峰值 明显 降低 ， 提 高 了 系统 的 响应 速度 ， 增 
强 了 系统 的 调节 能 力 ， 能 使 系统 很 快 达到 稳定 状态 ， 超 调 量 得 
到 了 降低 ， 提 高 了 系统 的 稳定 性 。 由 于 积分 阶 次 4 的 引入 ， 使 
控制 器 的 调节 更 有 具 灵 活性 。 优 化 后 得 到 的 控制 器 对 应 的 4 
=0.58， 低 于 图 像 法 得 到 的 4 =0.72. 由 于 4 的 降低 可 以 减 小 系统 
的 调节 时 间 ， 响 应 速度 快 。 通 过 仿真 验证 了 4 的 作用 ， 同 时 ， 
可 以 从 仿真 结果 看 出 优化 后 的 指标 明显 优 于 图 像 法 得 到 的 指标 。 

为 了 验证 优化 的 结果 为 最 优点 ， 分 别 在 w =16.43， 14 
=0.58 周围 进行 取 值 ， 因 为 其 附近 值 的 动态 特性 优 于 其 较 远 处 
的 值 。 所 以 仅 比 较 其 附近 的 值 。 值 的 选取 方式 如 下 : 

X1: w=14, 4=0.58 

X2: w=18, 4=0.58 

X3: w=16.43, 4=0.4 

X4: w=16.43, 4=0.7 
通过 第 2 节 中 的 公式 的 计算 ， 可 以 得 到 它们 对 应 的 控制 


C,(s) =2.76(1+ 3 


058 


CiG) = 2.88( + BP 


C;(s) =9.78(+ 路 


CuG) =5.03(1 + 7) 


然后 通过 式 (3)， 得 到 四 个 控制 器 对 应 的 系统 开 环 传递 函 
数 G3,G4,G5,G6. 则 它们 对 应 的 闭环 传递 函数 阶 跃 响应 曲线 如 
图 7 所 示 。 


Ampl 本 de 


图 7 不 同 控制 器 比较 
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通过 图 
标 ， 如 表 2 所 示 。 


形 可 以 得 出 不 同 控制 器 下 ， 系 统 阶 跃 响应 的 性 能 指 


不 同 控制 器 的 性 能 对 比 


调节 时 间 t， 上升 时 间 t， ” 超 调 量 o% 


G2 0.36 0.07 16.6% 
G3 0.44 0.078 17% 
G4 0.44 0.08 17.4% 
G5 达 不 到 给 定 值 
G6 0.5 0.09 18% 
右 下 角 是 图 形 ， 蓝 色 的 曲线 代表 C,(s) 用 人 工 鱼 


群 算法 优化 得 到 的 控制 器 ), 红色 的 曲线 代表 C,(s) ， 黄 色 的 1 


线 为 C,(s)， 紫 色 的 曲线 代表 C,(s) ， 绿 色 的 曲线 代表 C,(s) 。 


通过 表 2 和 图 7 看 出 ， C,(s) 的 调节 时 间 和 峰值 都 优 于 其 他 三 


个 控制 器 。 尽 


各 项 指标 ， 


为 了 验证 


兽 益 鲁 棒 1 


动 很 小 ， 但 达 不 到 稳 态 值 ， 综 合 考 虑 


所 器 最 优 。 


生 ， 分 别 取 ,=3.1383， ,=3. 8357， 


k, =3. 487 (分 别 是 标准 值 3. 487 的 +10%)， 通 过 取 不 同 的 k,， 


来 得 到 对 应 的 动态 响应 


Amplitude 


右 下 角 


线 ， 如 图 8 所 示 。 


二 Kp=3.1383 


一 一 Kpz=3487 


Kp=3.8357 


8 增益 鲁 棒 性 验证 


Time [s] 


红色 | 


以 看 出 他 们 的 民 
变化 为 +10%1 


则 满足 增益 鲁 棒 中 


4 


局 部 放大 图 形 ， 蓝 色 的 曲线 为 k, =3.487， 
线 为 k, =3.1383， 黄 色 曲 线 为 k, =3.8357。 通 过 图 像 可 


像 几 乎 是 重合 的 。 所 以 可 以 得 出 结论 ， 当 增益 
对 ， 超 调 量 几乎 不 变 ， 即 相位 裕 度 也 几乎 一 样 ， 
的 条 件 。 


二 阶 系统 分 数 阶 控制 器 的 整定 


对 于 二 阶 系统 ， 本 文 以 典型 的 二 阶 伺 服 系统 为 被 控 对 象 。 


近似 的 原理 


到 4. 则 此 时 P(S) 为 二 阶 传递 函数 。 对 应 


的 传递 函数 如 式 (22〉 所 示 。 


1 


+ (22) 


P(s) 


根据 系统 鲁 棒 性 和 稳定 性 的 要 求 ， 由 式 (4)~(6) 得 


[GGo,)|= 


k, Ja + ki @, cos(A5D)) + ki, sin(4757)) 


“(Yl + (OT ) 


(23) 


ChinaXiv 合 作 期 刊 


张 学 典 ， 等 : 基于 人 工 鱼 群 算法 的 分 数 阶 控制 器 参数 整定 


站 | _ kAw, sin(25) 
强 1 a 2 
dw vo ©. +2k@. cos()+k, (24) 
= 
1+(To.) 
kw” Sin(22) 要 
A = 1 Ew 2 1 T)= 
el 1 +kw.” oot I (25) 
-XxX/2+9, 
cot[tan (wo T) + 9,] 26) 


i 


@. “sin(2%) — @. “ cos(%)cotltan '(@, T) + 9,] 

给 定 截止 频率 、 相 角 裕 度 、 时 间 常 数 @ =5rad/s， mw, = 
60- ，7=0.5s。 分 别 画 出 式 (26)(12) 的 图 形 ，k 关于 4 的 函数 区 
像 。 如 图 9 所 示 。 


(26) 
(12)| J 


X: 0.593 
Y: 0.2425 


0 02 04 06 08 1 12 14 16 18 2 
lamda 


图 9 图 像 法 求 二 阶 系统 


从 图 像 可 以 看 出 (4 ,大 ) 的 值 为 《0.593,0.2425)， 代 入 式 
(23) 可 以 得 到 上, =2. 5484 进而 得 到 分 数 阶 Pr 控制 器 : 


C7(s) = 2.5484(1+ 和 (27) 


同时 根据 式 (21) 确定 寻 优 的 范围 ， 则 此 时 4 和 w 的 寻 优 
范围 为 0< wm, <10,0< 1 <1.186。 

根据 第 2 节 ， 用 人 工 鱼 群 算法 进行 寻 优 ， 则 可 以 得 到 w 
=6.04，4=0.41。 进 而 求 得 k; =0.0335, ,=3. 1314。 则 分 数 阶 
空 制 器 PT 控制 器 : 


谍 


Sa @8) 


分 别 求 出 对 应 的 二 阶 系统 的 开 环 传递 函数 G(s) ， G(s) ， 
对 它们 对 应 的 闭环 传递 函数 进行 仿真 ， 如 图 10 所 示 。 


Step response 


1.4 X: 1.35 
Y:1.265 a a 
‘2 NX: 1.49 G7 
X:0.91 312 X: 5.3 
人 1.004 Y: 1.016 
-加 
加 
包 
号 
豆 
三 
芯 
2 3 4 5 6 7 8 9 10 
Time [s] 


图 10 ”两 种 方法 的 对 比 图 
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把 两 个 系统 的 阶 跃 响 应 的 性 能 指标 以 表格 的 形式 进行 对 比 


Step response 


分 析 ， 如 表 3 所 示 。 era = 
表 3 两 种 方法 指标 对 比 ee | 
调节 时 间 太 上升 时 间 蕊 超 调 量 o% | s 
G7 5.3 0.91 31.2% Eb 
G8 4.3 0.86 26.5% Bo 
优化 得 到 的 分 数 阶 pp 控制 器 ， 调 节 时 间 ， 上 升 时 间 和 峰 %4[ 
值 降低 ， 提 高 了 系统 的 响应 速度 ， 增 强 了 系统 的 调节 能 力 ， 能 9 
使 系统 很 快 达到 稳定 状态 ， 超 调 量 得 到 了 降低 ， 提 高 了 系统 的 oo 0 
Time [s] 
稳定 性 。 
i 图 12 和 鲁 棒 性 的 验证 
同 理 在 其 周围 进行 取 值 ， 观 察 图 形 看 是 否 为 最 优点 。 
X5: w=5.04, 24.=0.41 通过 仿真 结果 可 以 看 出 : 不 同 的 kp 对 应 的 超 调 量 分 别 为 
X6: 四 =7.04，4=0.41 26.3%，27.6%，23.9%。 当 kp 变化 时 ， 系 统 对 应 的 超 调 小 幅 
X7: w=6.04, 4=0.5 度 的 变化 ， 对 应 的 相位 裕 度 也 相对 变化 较 小 ， 基 本 满足 鲁 棒 性 
X8: w=6.04, 4.=0.3 的 要 求 。 
通过 式 (23)~(27) 可 以 求 出 对 应 的 控制 器 : 人 
5 ”结束 语 


C,(s) = 2.5339(1 + 2) 


.0.41 
5 


首先 用 图 像 法 得 到 分 数 阶 pm 控制 器 的 参数 ， 在 此 基础 上 
指定 截止 频率 w 和 分 数 阶 系数 4 的 范围 ， 用 人 工 鱼 群 算法 对 
分 数 阶 PT 控制 器 参数 进行 优化 , 从 而 得 到 截止 频率 w 和 分 数 
阶 系数 1 ， 为 了 验证 用 人 工 鱼 群 算得 到 控制 器 最 优 ， 分 别 与 区 
像 法 得 到 的 控制 器 进行 仿真 对 比 ， 并 且 和 其 周围 的 值 形成 的 控 
制 器 进行 对 比 验 证 。 通 过 仿真 验证 可 以 看 出 ， 用 人 工 鱼 群 算法 
得 到 的 控制 器 动态 特性 最 优 。 为 了 验证 其 增益 鲁 棒 性 ， 通 过 对 

则 它们 分 别 对 应 的 二 阶 系统 的 开 环 传 函 数 。 比 k, 在 +10% 的 范围 变化 时 的 响应 曲线 ， 对 于 一 阶 系统 来 说 ， 
G9.G10.G11G12， 对 它们 对 应 的 闭环 传递 函数 进行 仿真 ， 如 “曲线 没有 明显 的 变化 ， 满 足 增益 鲁 棒 性 的 条 件 ， 而 二 阶 系统 的 
图 11 所 示 。 鲁 棒 性 没有 一 阶 系统 的 鲁 棒 性 好 ， 从 总 体 的 结果 来 看 ， 人 工 鱼 
群 算法 和 图 像 法 的 结合 使 得 系统 的 动态 性 能 指标 得 到 了 提高 。 


十 


CoG)=3.7804 -27) 


C1(s) =3.15( + 2%:) 


o 


CoG) =3.143( + 9% 


器 
y 


Step response 
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